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Abstract: 

This paper presents spatio-temporally resolved wideband measurements of Sub-Terahertz (Sub-THz) reflection coefficients in 
the frequency range of 92-110 GHz. A stochastic model for single reflection fixed links that is capable of modelling random 
scattering from small-scale discontinuities such as those encountered in complex structures in walls and partitions of buildings is 
presented. The model auto-regressively produces filter coefficients that are fed into an Infinite-Impulse-Response (IIR) filter which 
convolves them with the spatio-temporal series in order to generate the next output sample based on previous observations. 
The IIR filter allows for flexible stochastic generation of samples, and its parameters can be adjusted as needed to suit different 
channel conditions. A total of 20 suitable start-up filter coefficients are generated from 21.7 % of the sample size for each complex 
delay tap distribution, which corresponds to 50 complex instances of the channel. These coefficients are then utilised to validate 
the remaining measured sample set. The model is in quantitative agreement with measurement statistics and can be used to 
construct relatively simple modified reflection coefficients that can be used in micro-cellular ray-optical network planning tools. 


1 Introduction where the direction of rays are suitably defined [8]. Similarly, util- 


Recently, there has been significant interest in the Sub-Terahertz 
(Sub-THz) band, which covers frequencies ranging from 0.1 to 
0.3 THz. With its wide spectrum availability, communication sys- 
tems operating in the Sub-THz band can support more users and 
deliver faster data rates [1]. This frequency range has attracted 
attention due to its potential for various applications, including high- 
capacity fixed links in indoor locations to overcome difficulties 
with employing fibre [2]. This can be primarily accomplished by 
providing Line-of-Sight (LoS) fixed links that are comparable to 
Free Space Propagation (FSP) [3]. However, establishing directed 
Non-Line-of-Sight (NLoS) fixed links with a single reflection is 
crucial for enabling effective communication and overcoming the 
limitations imposed by the high path loss and directivity of the 
antenna. With the advancements in digital wideband communica- 
tions, it has become imperative to incorporate the spatio-temporal 
and frequency selective characteristics for an accurate representa- 
tion of scattered communication channels for both mobile and fixed 
radio communications [4]. 

Micro-cellular radio planning tools at Sub-THz frequencies are 
primarily based on Ray-Tracing (RT) algorithms [5]. Propagation 
studies in micro-cellular environments have shown that significant 
multipath components arise from reflections of building surfaces [6]. 
Reflections from indoor building surfaces at Sub-THz are frequently 
modeled as “homogenised” lossy dielectric planar surfaces. Such an 
assumption permits the use of simple Fresnel reflection coefficients 
to be employed in modelling the reflection of electromagnetic waves 
in a computationally efficient manner. This, however, is only true 
for the specular reflected path. In reality, scattered wave distribu- 
tions resulting from propagation through ’smooth’ looking complex 
structures that contain features such as metallic reinforcement studs 
and insulating material are highly complex, as they vary significantly 
from one point to another within the structure [7]. Attempting to 
model such complex wave distributions solely based on geometric 
considerations becomes inadequate. As a result, such distributions 
cannot be modeled as Rayleigh or Rice as per e.g., Jakes’ model, 
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ising Kirchhoff scattering [9] theory to model these distributions 
based solely on the surface roughness proves insufficient, as it over- 
looks the internal scattering effects. While existing literature has 
extensively characterised homogeneous [10] and stratified [11] sur- 
faces, there is a notable gap when it comes to the characterisation 
of surfaces with random internal structures at Sub-THz frequencies. 
In some cases, these surfaces may be smooth on the outside but 
cause substantial multipath at these wavelengths due to their inter- 
nal structure. Conventional RT methods are not sufficient to model 
such characteristics. Hence, a new approach is required to model the 
spatially-variant and frequency selective characteristics of randomly 
scattered communication channels based on the displacement along 
the surface. 

To generate spatio-temporal fading estimates for statistically cor- 
related delay taps, it is assumed that the propagation channels 
enabled by Single-Input-Single-Output (SISO) propagation with 
very localised scattering are of Wide-Sense-Stationarity (WSS) [12]. 
Under this assumption, the fading on each delay tap is treated as 
independent and is modeled as a displacement-series process. How- 
ever, collecting the necessary data to capture the underlying propa- 
gation characteristics can be challenging, as it requires a significant 
amount of time and effort. Therefore, due to practical limitations, 
only a limited set of data can be acquired during the measurement 
process. Given these challenges, recent advances in machine learn- 
ing that encompass techniques such as GANs (Generative Adver- 
sarial Networks) [13] which allow for complex distributions to be 
captured, may be difficult to train and computationally inefficient. 

In this context, we present experimental spatio-temporally 
resolved measurements of wideband continuous channel coefficients 
at 92-110 GHz, to quantify scattering due to small-scale discontinu- 
ities found in complex structures. We explore the application of an 
Auto-Regressive (AR) Infinite-Impulse-Response (IIR) filter-based 
model to formulate an equivalent frequency variant reflection coeffi- 
cient that can be used in RT tools. By estimating the AR coefficients 
and convolving them with the input data using an IIR filter, the model 
can generate predictions for future values of the displacement-series 


Fig. 1: Internal scattering from complex building structures for the case 
where (0; = 0r). 


based on its past values. To the best of our knowledge, no study 
on the application of the AR approach at Sub-THz frequencies to 
accurately quantify spatio-temporally resolved propagation through 
complex structures has been reported in the literature. 


1.1 Multipath from Complex Internal Building Structures 


At Sub-THz frequencies, the presence of complex building struc- 
tures gives rise to multipath propagation, which is primarily 
attributed to random scattering phenomena occurring within the 
structures themselves. At these frequencies the wavelength becomes 
shorter, resulting in a smaller distance between consecutive wave 
crests or troughs. When electromagnetic waves encounter small- 
scale discontinuities within the structures, such as irregularities or 
changes in materials, they undergo random reflection and dispersion 
in various directions. Consequently, a complex propagation pattern 
emerges inside the structure that cannot be accurately modeled using 
simple geometric principles. 

These complex building structures encompass various elements, 
such as plastered walls containing symmetrically placed vertical 
studs and insulating material positioned between the studs. Fig. 1 
provides a visualisation of the internal scattering effects when the 
incident angle (0;) is equal to the reflection angle (0r). Notably, the 
directed wave undergoes random reflections from different points 
within the structure, leading to delays in the arrival of multipath 
signals and variations in their power levels. Overall, the complex 
nature of Sub-THz multipath propagation arising from random scat- 
tering within complex building structures necessitates a more sophis- 
ticated modelling approach that goes beyond external geometric 
considerations. 


1.2 Structure of the paper 


The remainder of this paper is organised as follows; Section 2 details 
the measurement setup, calibration, and data collection processes. 
Section 3 involves the analysis and appropriate modelling of the col- 
lected data, which is subsequently utilised in the training algorithm. 
The model for generating fading on individual delay taps is then 
formulated in section 4 and its performance is verified against the 
measurement data. Section 5 concludes the findings presented in this 


paper. 


2 Measurement Setup & Calibration 


The measurement campaign setup is shown in Fig. 2. A N5230A 
Vector Network Analyzer (VNA) from Keysight was used to obtain 
the S21 scattering parameter, operating in the 2 to 20 GHz frequency 
range and up/down converted. The Tx and Rx sides consist of a Viva- 
Tech THz up-converter and down-converter module, respectively, 
to achieve the desired frequency range of 92-110 GHz. The fre- 
quency extender modules are then connected to WR-10 Waveguide 
test ports, and two nominally identical VDI WR-10 diagonal horn 
antennas of 26 dBi gain and a Half-Power-Beam-Width (HPBW) of 


7.5° at the center frequency, one for each of the extender heads at 
the Transmitter (Tx) and Receiver (Rx) ends, to transmit the Sub- 
THz signal. The scenarios under investigation seen in Fig. 2 include 
scanning over a plastered wall with insulating material behind (a), 
scanning over a metal stud found behind the plastered wall (b), scan- 
ning over the plastered wall while approaching a metal stud (c), and 
scanning over a wooden panel (d). The studs for measurements at 
locations (b) and (c) are type °C?’ studs with a known width of 48 
mm that were previously identified using a metal detector. The Tx 
and Rx were placed at the edge of a platform controlled by a linear 
mechanical positioner sitting on a trolley, at a 1 m height and a 10° 
specular configuration where the total propagation distance from the 
Tx to the reflector to the Rx was 2 m. The linear measurements were 
calibrated beforehand by placing a metal sheet in front of the reflect- 
ing surface to minimise any losses due to measurement equipment. 
A total of 201 data points were acquired during frequency scanning 
for each location. Given an 18 GHz bandwidth, this corresponds to 
a distance resolution of 16.7 mm between consecutive delay taps 
in free space, and a maximum excess delay of 3.35 m, sufficient to 
capture any multipaths originating from the region of interest. 

For each scenario, a total of 230 frequency scans were acquired 
during the measurement using a linear mechanical positioner with a 
resolution capability of 1 mm. Provided a center frequency of 101 
GHz, this calculates to give a wavelength of ~ 3 mm. The setup was 
shifted by As, which corresponds to a change in specular position 
on the reflective surface. As should be small enough that it does not 
move beyond a portion of a wavelength from sample to sample, so 
that the channel characteristics do not vary significantly from point- 
to-point. Therefore, the sampling resolution must be in excess of 
twice the Doppler shift of the channel. By forming a Doppler deriva- 
tion, an appropriate sampling rate for the measurement campaign 
can be determined. The Doppler period, Ty can be described as the 
ratio between the wavelength, A at the centre frequency and velocity 
of the specular position, Vs given by: 


— Àc 
VW 
As the Doppler derivation must be established in a way were the 
samples are taken over displacement instead of time since our mea- 
surement is static and there isn’t a corresponding velocity, samples 
can be taken relative to the wavelength using factor k, as follows: 


Ta a) 


eee 2 
wae (2) 
The sampling period can then be formed as: 
Ac As 
T. = ze 3 
si, = V (3) 


which can then be normalised for any given velocity such that then 
the following equation must hold: 


—1_& As 
fs k Ta Ac 


Therefore as it is necessary to increase the sampling frequency rela- 
tive to the Doppler frequency by factor k, where 2 is a suitable choice 
according to the Nyquist-Shannon sampling theorem, then samples 
must be taken every half wavelength [11]. Thus, As must not be 
greater than 1.5 mm at every sample point. In this case, we select 
factor k to be equal to 3 and sample every 1 mm to achieve a high 
measurement auto-correlation in order to enhance the capabilities of 
the AR process. 


(4) 


3 Spatio-Temporal Measurements & Analysis 


This section presents the time-domain analysis of the measured data 
for the scenarios considered. The Power Delay Profiles (PDPs) are 
initially presented and the delay taps of interest are further analysed 
to evaluate the effects of propagating through the structures under 
investigation. 
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Fig. 2: Spatio-temporal reflection coefficient measurement setup. (a) Between wall studs. (b) Over metal stud. (c) Approaching metal stud. (d) Over a 
wooden panel. 
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Fig. 3: Power delay profiles. (a) Between wall studs. (b) Over metal stud. (c) Approaching metal stud. (d) Over a wooden panel. 
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Fig. 4: Multipath power distributions over scanning range. (a) Between wall studs. (b) Over metal stud. (c) Approaching metal stud. (d) Over 
a wooden panel. 


3.1 Power Delay Profiles specular point on the surface and arriving at delay To = 6.67 ns, cal- 

culated given the speed of light in free space and a total propagation 
Sample PDPs for the four scenarios considered are shown in Fig. 3. distance of 2 m. For all scenarios, the PDP has been filtered so that 
The PDPs were normalised so that the the strongest component at it includes delay distributions with spatial components that at any 
each spatial offset point is offset to 0 dB magnitude, then all other point during measurement reached above -30 dB while values below 
samples are offset in magnitude by the same value. Usually, the are set to -70 dB for clarity. Fig. 3 (a), depicts the case where fre- 
strongest component for all scenarios is the one reflecting from the quency scans were taken over a wall between vertical metallic studs 
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placed symmetrically along the structure. Between these studs, insu- 
lating material can be found with high non-homogeneity and a lot 
of random scattering can be expected from various points internally. 
By first glance, randomly scattered components are observed arriv- 
ing after the specular peak of the PDP that reach up to about -10 
dB below the peak. Fig. 3 (b), illustrates the cases were frequency 
scans were taken over a metal stud and beyond. As expected a sec- 
ondary strong delay distribution can be observed with components 
reaching the specular components in magnitude over the range of 
40-88 mm and then dropping, which verifies the presence of the 
metallic stud internally within the structure. After scanning over the 
stud, a second peak seems to appear before the magnitude further 
drops. This may be due to multipath being picked off the edge of 
the stud. Fig. 3 (c), illustrates the case were frequency scans were 
taken over the complex wall structure while approaching a metal 
stud. Strong components can be observed for the first 100 mm due 
to reflections occuring internally that may be due to wires or sock- 
ets, while the magnitude after about 210 mm reaches the specular 
peak as in case (b), indicating that another stud has been reached. 
Fig. 3 (d), demonstrates the case where frequency scans were taken 
over a wooden surface. Due to the fact that this type of surface is 
considered to be highly homogeneous, much less randomness in 
scattering is expected along with a flatter and less variant response 
in the frequency domain, unlike the measurements conducted over 
the complex plastered wall structure. 

By further inspection of multipath distributions over the scanning 
range as shown in Fig. 4, significant distributions are identified and 
analysed in relation to the structure of the surfaces investigated. Any 
multipath components below -20 dB have insignificant contributions 
to the overall power and can be treated as noise. For simplification 
purposes, the first two tap distributions that arrive after the reflec- 
tion of the plastered wall and before the stud reflection for cases 
(b) and (c) at To + 0.55 ns and 79 + 0.11 ns which are mostly in the 
noise, have been omitted for simplification purposes. Similarly, these 
taps have also been omitted for cases (a) and (d) as they also appear 
weak in magnitude and therefore the model statistics can be com- 
pared fairly against all scenarios by considering the same number 
of delay taps. In Fig. 4 (a), the distribution of multipaths over scan- 
ning distance are observed for the between studs case. The taps that 
follow the specular peak include strong components arriving from 
the structure internally that reach above -20 dB at 79 + 0.165 ns and 
To + 0.22 ns, followed by some weaker more complex distributions 
that fluctuate heavily between -20 to -40 dB. In Fig. 4 (b), the dis- 
tribution of multipaths can be observed for the case were frequency 
scans were taken over and beyond the stud. The stud region is high- 
lighted in the figure for clarity. The delay component that comes off 
the stud is seen to arrive 0.165 ns after the main peak. This delay 
corresponds to an additional propagation distance of approximately 
33.34 mm, by taking into account the measured permittivity of the 
plasterboard, er = 2.2 [14] at 100 GHz and the speed of the wave 
propagating through the material [15], which is calculated as: 


c 

w= (5) 
This suggests that the plasterboard thickness is about 16.7 mm, while 
the resolution of the delay tap is every 11 mm in the plasterboard so 
it is reasonable that it is detecting the stud even though there is some 
uncertainty of this value. Furthermore, despite this uncertainty, the 
delay components following the one that relates to the stud at 7) + 
0.22 ns, To + 0.275 ns and To + 0.33 ns, though lower in magni- 
tude, still pick up off the stud as the distributions look similar to that 
of the stud at 79 + 0.165 ns. This is probably due to multiple reflec- 
tions occuring within the air interface between the stud and plastered 
wall as the additional delays beyond the stud correspond to propa- 
gation distances greater than the thickness of the plasterboard. Fig. 
4 (c), illustrated the distribution of multipaths while approaching a 
metal stud at 210 mm where the stud component at To + 0.165 ns 
approaches the strength of the specular component. This indicates 
that more energy is passing through the plasterboard and stronger 
reflections are arriving from the insulating material. Interestingly, 
before the stud and between about 50-100 mm multiple reflections 
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Fig. 5: PDP correction factor relative to the specular component 7o. (a) 
Between wall studs. (b) Over metal stud. (c) Approaching metal stud. (d) 


Over a wooden panel. 


appear internally due to the distributions at 7) + 0.22 ns and To + 
0.275 ns looking very similar to that of the stud picking off the at 
To + 0.165 ns. These reflections may come from wires and sockets 
internally as they are consistent and are possibly not coming off the 
insulation material. Lastly, in Fig. 4 (d) multipath distributions can 
be observed when scanning over a wooden panel. Besides the first 
path arriving at 79, The taps that arrive at To + 0.165, To + 0.22 and 
To + 0.275 are significant since they reach well above -20 dB while 
the remaining two taps that are much more randomly scattered and 
are below -20 dB can be treated as having insignificant contributions. 

To recover the un-normalised PDP, a correction factor can be 
applied as shown in Fig. 5, which essentially illustrates the specu- 
lar component, 79, before the normalisation process. Importantly, a 
-6 dB drop can be observed where there is a stud for measurement 
cases (b) and (c). This is because the stud only allows two paths to 
go through, but if more paths are added in addition to the original 
specular path from the plasterboard when away from the stud, then 
that means the mean magnitude of the channel has to increase if taps 
are added. Thus, this is seen when away from the stud. So, it can be 
concluded that the stud locations are now known. This means that 
the model can accommodate this correction factor, which for studs 
can be an adapted periodic function if the studs are positioned with 
equal spacing. For the other cases, (a) and (d), the correction factor is 
not as significant as the delay taps are more consistent in magnitude. 


4  Auto-Regressive IIR Filter based Model of 
Small-Scale Fading 


This section presents the application of an AR-IIR filter based model 
to characterise the scattering observed in the PDPs presented in 
section 3. The model is further validated against the first and sec- 
ond order statistics of the propagation channel and the channels are 
finally reproduced in the frequency domain, aiming to capture the 
frequency selective effects of such displacement-variant channels. 
AR models are commonly used for time series prediction, and IIR 
filters can be applied to implement an AR model. The AR process 
employs information from preceding states as input to a regres- 
sion equation to anticipate the sample at the subsequent time step. 
AR models have previously been employed for Kalman-filter-based 
channel estimation [16], [17] and for channel propagation predic- 
tion in the lower sub-6 GHz frequency bands [18], [19]. They have 
also been thoroughly investigated for the synthesis of correlated 
Rayleigh fading bandlimited processes [20]. AR modelling depends 
upon a limited number of parameters such as the ones acquired 
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Fig. 6: Framework of the proposed auto-regressive IIR filtering model. 


during measurement. A number of different methods exist to pre- 
dict the AR parameters such as the Yule-Walker, least-squares and 
Burg’s method [21]. The least-squares and Yule-Walker methods, 
which estimate the AR parameters directly cannot be employed, 
as the model generated is unstable yielding poor parameter esti- 
mates [22]. Burg’s method on the other hand, initially provides an 
estimate of the reflection coefficients that give the last parameter pre- 
diction for each past observation p where P is the last observation in 
the displacement series. This is done by minimising estimates of the 
forward and backward prediction error powers for past observations 
in the least squares sense, recursively, and are constrained to satisfy 
the Levinson-Durbin recursion algorithm. When the reflection coef- 
ficients generated based on the input displacement series become 
small in magnitude compared to the preceding ones, this indicates 
an appropriate model order for the AR process. 

The framework for the proposed AR-IIR filtering model is pre- 
sented in Fig. 6. The AR process starts by initialising the estimated 
coefficients to zero and then updates them one by one. At each iter- 
ation, the algorithm uses all the available data points to estimate the 
next coefficient. This is due to the fact that the algorithm needs a 
sufficient number of data points to ensure accurate estimation of 
each filter coefficient, that is not necessarily equal to the number 
of generated filter coefficients, suggesting a trade off between accu- 
racy and complexity. Based on estimates of the reflection coefficients 
Ky, alil, Kr;2[9], .-- K7, p[j], the AR filter parameter estimates 
Gr, 1l], Gr;,2[9], ---. âr p[i], and prediction errors ê+; [j], for 
each delay tap 7;, and for j = 1, 2,..., J iterations signified by the 
number of required predictions, are estimated recursively as follows: 


Forj = 1: 


ân pll] = 1 for p = 0 
TPES T Š- pll] for p = 1,2,..., P 


ên [1] = oh, - (1 = lân; Pill?) 


where oF, represents the variance of the input sequence. Based 
on the initial filter coefficients @7,,»[1] and initial prediction error 
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é7,[1], the next set of coefficients G7; ,»[j], and prediction errors 
êr; [j], are generated iteratively as follows: 


For j = 2,3,..., J: 


1 for p=0 
Gr: pli] = § âri plil + Kr: pli] 47,919 - 1] forp=1,...,P-1 
Kz plj] for p = P 
(6) 
ên [j] = ên [j — 1] - (1 — lân, Plil’) (7) 


The first coefficient G;, 0, is always set to 1, which is a common 
initialisation practice in AR modelling. The last filter coefficient 
prediction @,, p equals to the last reflection coefficient prediction 
an p for the je iteration as there are no more filter coefficients 
to update. The filter coefficients for p = 1 to P - 1 are generated 
based on previously computed filter coefficients, âr; ,p[j - 1] to obtain 
the current set of filter coefficients, â+; ,p[j]. Similarly, this recursive 
process uses the previously computed prediction error êr; [j — 1] to 
obtain the current prediction error ê+, [j]. The initial values of the fil- 
ter coefficients â; [1] and the prediction error, ê+; [1] are therefore 
updated throughout the iterations to refine the estimates. 

The IIR modelling process builds on the AR parameters esti- 
mated through Burg’s method to generate the next successive sample 
predictions as follows: 


he, pj =- >, Gri plil * n, P-ptj + érfi]) 8) 


The equation represents the recursive prediction process of the out- 
put of the Infinite Impulse Response (IIR) filter, denoted as h+, p+, 
for each delay tap 7 and for J predictions, where P in this case not to 
be confused with the total number of filter coefficients, is used an as 
index that indicates the last observation. The prediction at each step 
j depends on the AR filter coefficients @7; ,p and the prediction error 
Pr;,p Obtained at model order P. The recursion involves convolving 
p = 1 to P filter coefficients and the previous observations including 
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Fig. 7: Validation of delay tap distribution at To + 0.275 ns. 


the predicted ones hr, P—p+ lJ — 1] summed up with the predic- 
tion error €7;,p in each iteration to ensure that the IIR filter takes 
into account an equal number of filter coefficients and past observa- 
tions to make the predictions for each delay tap 7;. The negative sign 
in the equation signifies the feedback nature of the IIR filter, where 
the outputs are subtracted from the current input, creating an infi- 
nite impulse response. The IIR modelling process provides valuable 
insights into the time-domain behavior of the channel’s multipath 
components, contributing to a more accurate representation of the 
underlying propagation channel. 

The first few filter coefficients generated are strong in magnitude, 
which implies that the remaining weaker ones are used for fine tun- 
ing. A reasonable number of P = 20 filter coefficients generated from 
the first 50 complex channel instances are used for the training pro- 
cess to make the first prediction, that is then padded to the predicted 
series to generate the next set of coefficients. This process occurs 
iteratively, where the filter coefficients are generated based on all 
sample data points including an additional sample point that is pre- 
dicted by the IIR filter in each iteration, to eventually validate the 
remaining sample size that consists of 180 reflection coefficients. 
By generating tuned start-up filter coefficients that give a reason- 
able response, the stability of the model can be improved as an IIR 
filter may easily diverge or provide random results with improper 
coefficients. 


4.1 First Order Statistics 


In order to validate the proposed model and its suitability for prac- 
tical use, its performance is assessed by comparing it against the 
small-scale characteristics of the propagation channel using mea- 
sured data. These statistics are instrumental in characterising the 
overall distribution of the fading multipaths. For illustration pro- 
poses, one of the most complex and significant distributions is 
selected, specifically the one that arrives 0.275 ns after the specu- 
lar component, To. The small-scale attenuation of the propagation 
channel can be described in terms of first-order statistics, such as the 
magnitude of the channel impulse responses as seen in Fig. 7 which 
can be also used to derive the Cumulative-Distribution-Function 
(CDF) as illustrated in Fig. 8. In all cases there is good agreement 
between the measurement data and model without any significant 
discrepancies. Nonetheless, since the AR process can produce com- 
plex outputs, it also allows for the angular component to be resolved 
which is simply a phase gradient, and is therefore omitted from 
illustration. 
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Fig. 8: Validation of the cumulative distribution for the proposed frame- 
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Fig. 9: Validation of the Doppler spectrum at To + 0.275 ns. 


4.2 Second Order Statistics 


While the first-order statistics provide valuable insights into the over- 
all distribution of fading multipaths, they don’t convey information 
about how these multipaths are distributed over time. In addition to 
the first-order statistics, the second-order statistics are also derived, 
such as the Doppler spread, Level-Crossing-Rate (LCR) and Aver- 
age Fade Duration (AFD). These statistics provide information on 
how fading multipaths are distributed over time. This information 
is valuable for designing and evaluating communication systems. 
Fig. 9 illustrates the Doppler spread that verifies the suitability of 
the model for small-scale characterisation of each measurement sce- 
nario. The x-axis represents frequency relative to the maximum 
Doppler shift, fa vax, Which ranges from - 1.5 fa Max to+ 1.5 fa Max: 
The peak magnitudes show approximately zero Doppler shift, and 
the Doppler components fall to the noise floor more than 20 dB 
down from the peaks. Narrow peaks can be observed due to the 
small amount of measurement data, while any transients around the 
peaks are well captured before the noise floor is reached. Fig. 10 
validates the threshold for level crossing against the LCR, which 
can be expressed as the percentage of times per second the fading 
signal falls below a certain level threshold relative to the Doppler 
period (1/ fa). A further validation of the second order statistics is 
presented in Fig. 11, which shows a good agreement between level 
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Fig. 10: Validation of the second order statistics based on level crossing 
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Fig. 11: Validation of the second order statistics based on average fade 
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crossing rate and AFD which can be expressed as the average time 
the fading signal is below a certain level threshold relative to the 


Doppler period (1/ fa). 


4.3. Frequency Domain Analysis 


The effectiveness of the proposed approach in recovering the fre- 
quency selectivity of the channels is demonstrated through frequency 
domain analysis. In Fig. 12, plots of the complex channel coefficient 
versus frequency are presented for all investigated scenarios. The fig- 
ures depict the response at three specific prediction points: the first 
prediction point (offset point 51), the middle prediction point (offset 
point 140), and the end point of the prediction (offset point 230). The 
results clearly illustrate that the proposed approach accurately cap- 
tures frequency selective dips in power and effectively tracks shifts 
in frequency selectivity. The plastered wall structure investigated 
clearly forms much larger dips in frequency selectivity compared 
to the much more homogeneous response of the wooden panel, 
which is attributed to the weaker effects of scattering. These find- 
ings indicate the robustness and reliability of the proposed approach 
in recovering the frequency selective characteristics of the scattered 
communication channels. 


IET Research Journals, pp. 1-8 
© The Institution of Engineering and Technology 2015 


5 Conclusion 


In summary, the application of a simple Auto-Regressive IIR filter- 
based model for characterising small-scale discontinuities compara- 
ble to the wavelength resulting from reflections of complex struc- 
tures is proposed. An overview of the methodology outlining the 
correct measurement practices and procedures for spatio-temporally 
resolved scattered propagation was presented. This was followed by 
a Statistical comparison of the measurement results with the model 
predictions that were used to verify the applicability of the model 
for use in RT planning tools. Drawing upon the results presented, 
preliminary conclusions can be made regarding the implications for 
system design and deployment: 


e In the context of modelling propagation through complex struc- 
tures, RT tools can leverage the periodic nature of these structures. 
This involves individually modelling specific sections of the struc- 
tures, which enables the formation of a Markov chain of statistical 
distributions. 

e The presence of studs within complex structures give rise to an 
additional strong path beyond the specular reflection arriving off 
the stud, and further complex distributions coming off the internal 
structure that can be accurately captured by the AR process. 

e Considering the effectiveness of the proposed model, it is neces- 
sary to obtain tuned start-up filter coefficients to stabilise the model 
as an IIR filter may easily diverge or provide random results with 
improper coefficients. 


The drawn conclusions offer valuable insights into the planning 
and deployment of Sub-THz systems. The practices showcased in 
this paper and the resulting findings hold substantial interest for spe- 
cialists and standard bodies engaged in propagation modelling for 
the design and deployment of communication systems. 
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